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I. INTRODUCTION 



The growth of ocean acoustics as a theoretical and 
experimental science can be traced to a fundamental change 
in the conduct of warfare. During World War II, the 
belligerents employed scientific technology for the conduct 
of war to a degree’ theretofore unknown. In the United 
States scientific research was mobilized on a scale 
comparable to the mobilization of manpower and industry. 
Whole fields of scientific research were expanded to produce 
a technology for the conduct of the war. A portion of this 
research, the study of propagation of sound in the sea, was 
a response to a particularly effective (and technical) enemy 
weapon - the U-boat. 

The first wartime studies of underwater sound 
propagation hinted at sound transmission over long ranges 
within a natural oceanic waveguide known as the SOFAS 
channel. A series of experiments was conducted in which 
explosive signals were propagated in the SOFAS channel over 
distances of up to 900 nautical miles (nm) [ 1 ]. In 
subsequent experiments explosive signals were propagated 
over paths of up to 2,300 nm [2]. Since World War II 
several uses for SOFAS propagation have been devised, 
including the location of missile impacts within wide 
oceanic areas and the study of the low-frequency ambient 
noise field. 

To predict and aid in understanding SOFAS propagation, a 
number of theoretical models for the transmission of sound 
in the ocean have been developed. The theoretical 
treatments have fallen generally into two broad catagories: 
ray-tracing and normal mode techniques. Ray-tracing 
techniques are based on the Eikonal equation. 
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( 1 . 1 ) 



2 2 
V A • VA = N = [ C /C (Z) ] 

o 

where A is the ray function, N an index of refraction, c is 

o 

a reference sound speed, and c(z) sound speed as a function 
of depth. The Eikonal equation is restricted to high 
frequencies in which the sound speed fluctuation is small 
over one wavelength. Normal mode techniques, more suited to 
low frequency sound, are based upon the depth* dependent 
Helmholtz equation, 

2 2 2 

[ V + a) /c (z) ] y = 0, (1.2) 

where Y is the acoustic velocity potential and go is the 
angular frequency of the acoustic disturbance. 



A. DEVELOPMENT OF THE NORMAL MODE TECHNIQUE 

Since the work to be presented is based upon the normal 
mode technique, a brief overview of the development of 
underwater sound transmission models based upon normal mode 
theory is in order. 

1 • Early Theo ry 

During World War II Pekeris developed two 
mathematical approaches which became the foundation for much 
of the later work. In 1946 Pekeris published a treatment of 
sound propagation in a layer characterized by a variety of 
sound speed gradients [3]. Although he was primarily 
concerned with sound penetration into the shadow zone caused 
by a negative sound speed gradient, he devleoped asymptotic 
solutions to the wave equation for a number of sound speed 
gradients. For the case of a linear sound speed gradient, 

c = c (1 ♦ az) , (1.3) 

o 
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he derived a solution which reduced to a Fourier- Bessel 
Series, which converged rapidly in the shadow zone and 
slowly in the insonified zone. In addition, he investigated 
a number of variable sound velocity profiles which included, 

2 2 -1 

c = c (1 - az) , (1.4) 

o 

for which the solution is a combination of Hankel functions 
of order one-thitd. 

In 1948 Pekeris published a method of determining 
the normal modes and resultant sound field in both 

two-and-three-layer isovelocity media [4]. To evaluate the 
sound field he applied a Hankel transform in the complex 
horizontal wave number plane to a vertical wave function U, 
which resulted in an expression for the acoustic velocity 
potential 

Y = fuj k dk . (1.5) 

c r o r r 

The wave function displayed a number of singularities and a 

branch cut. The solution consisted of a summation of normal 

modes, represented by a series of simple poles, and a ground 

wave in the lowest layer, represented by a branch cut (see 

Fig 1) . The branch cut extended from the branch point 

corresponding to c , the sound speed in the bottom, to i 

b 

This particular cut necessitated the evaluation of 
additional complex- valued poles between the branch cut and 
the positive imaginary axis (the so called "leaky modes") . 
The "leaky" modes were seen to be important only at short 
ranges. Pekeris evaluated the integral by the calculus of 
residues for the poles and the evaluation of a branch cut 
associated with the sound speed in the lowest layer. 

There were two significant modifications to the 
Pekeris technigue which soon followed. Tolstoy reformulated 
the characteristic equation for the modal eigenvalues 
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FIGUBE 1 - Path of Bessel Integrals in Complex Plane 



16 



(location of the poles in the complex wave number plane) in 

terms of reflection coefficients for both upgoing and 

downgoing waves [5]. Tolstoy considered the effects of 

horizontally stratified layers above and below a depth of 

interest, z . A wave travelling up past z will be 
o o 

reflected with a reflection coefficient Bf. Similarly a 

plane wave travelling downward past z will have a 

o 

reflection coefficient R 4 . . Tolstoy showed that for an 

undamped mode representing the propagation of energy in the 
waveguide, the reflection coefficient above and below z 

o 

must be related by 



Bf R* = 1 . (1.6) 

Rf and Rf are functions of the sound speed profile and a 
horizontal parameter (such as the horizontal wave number or 
the horizontal phase speed). Since Eg. (1.6) is satisfied 
only for certain discrete values of this parameter, it 
represents a characteristic equation for the modal 
eigenvalues. This general characteristic equation enabled 
normal mode analysis to be applied to a wide variety of 
environmental conditions. 

The similarity between the normal mode problem in 
acoustics and the SchrOdinger equation in quantum mechanics 
was recognized and provided a tool for further development. 
Bremmer published an analysis of the Wentzel-Kramer- 
Brillioun (WKB) method of quantum mechamics as an 
approximation to wave propagation in an inhomogeneous 
medium, opening the way for the use of the WKB solutions and 
characteristic equations as solutions to the vertical wave 
equations [ 6 ]. 

The second change to the Pekeris technique came when 
Ewing, Jardetzky and Press [7] suggested an alternate path 
for the branch cut (Fig 1). The Ewing, Jardetzky and Press 
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branch cut extended from the same branch point, to the 
origin, and then along the imaginary axis. This eliminated 
the need for evaluation of the complex-valued poles, as they 
were shown to be displaced to an adjacent Riemann sheet. 

The Pekeris model and its descendents were used in 
many theoretical applications, including shallow water 
propagation over fluid and elastic bottoms [8 - 15]; and was 
also used to study the effects of attenuating [16], sloping 
[17], and rough bottoms [18] on mode propagation. 

Up to this point a major limitation for normal mode 
and many ray-tracing techniques was the required assumption 
that sound speed be a function onl£ of depth (to separate 
the range and depth-dependent equations) . In the ocean the 
sound speed profile varies gradually (and at times abruptly) 
with range. The applicability of the normal mode technique 
and, to a lesser extent, the ray-tracing technique over a 
horizontally varying medium was broadened by Weston in 1958 
with the development of mode and ray invariants [19]. 
Weston observed that a mode may be assumed to keep its 
identity for slow horizontal changes in the sound speed 
profile. This is known as the adiabatic assumption (after 
its analogue in time-dependent perturbation theory of 
quantum mechanics) , implying that the energy propagated by a 
particular mode remains associated with that mode. A more 
receat paper by Milder expands upon the concepts of mode and 
ray invariants and derives a test for the applicability of 
the adiabatic assumption [20], The adiabatic assumption 
allows local horizontal stratification in an environment 
which changes slowly in the horizontal. Thus the normal 
mode method now became applicable to a wider variety of 
situations provided the horizontal oceanic variation was 
gradual over the long ranges. In addition, although not 
strictly required, the assumption of local horizontal 
stratification has been found useful in some ray-tracing 
work [21]. 
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2. Num e rica l Meth ods 



In the late 1950's and early 1960's there began to 
appear a number of numerical normal mode programs which took 
advantage of the digital computer. Many of the models were 
based at least in part on known solutions for certain 
classes of sound speed profiles. 

In 1959 Tolstoy and May developed a normal mode 
computer model based upon a series of horizontal layers in 
which the inverse of the sguare of the sound speed varied 
linearly with depth 



c = c (1-az) . (1.7) 

o 

In such layers Tolstoy and May approximated the mode 
solutions (Bessel functions of order one-third) with 
twelfth-order polynomials [22]. For a given horizontal wave 
number they started with both a surface and bottom boundary 
condition, then iterated from the boundaries to some 
interior depth. The Bessel functions were matched between 
layers by invoking continuity of acoustic pressure and 
vertical particle velocity. At the interior depth the 

slopes of the solutions iterated from the surface and from 
the bottom were compared. Only certain values of the 
horizontal wave number, which defined the eigenvalues, gave 
first order continuity across the interior depth. The 
solutions thus formed defined the normal modes. 

Peter Hirsch studied axial sound focusing in the 
SOFAB channel using an analytic representation for the sound 
speed profile [23]. He assumed a sound speed profile 
-2 

parabolic in c and symetric about the channel axis, 

2 2 22-1 

c = c (1-a z ) . (1.8) 

o 

For a pulse source this representation of the SOFAB channel 
gave successive focusing at about 6 nm intervals. Williams 
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and Horne did another study of SOFAR channel axial focusing 
using a polynomial representation for the sound speed 
profile [ 24 ], 

2 2 2 3 4 -1 

c = c (1-az-bz -cz -dz ) . (1.9) 

o 

This polynomial representation of the sound speed profile 
gave a better fit to the actual SOFAS channel, allowing 
representation of its asymmetrical character. Solving for 
the eigenvalues and normal modes, they described the sound 
field for a continuous wave (CW) source. Williams and Horne 
found a much smaller degree of axial focusing than Hirsch; 
they also determined that focusing effects were sensitive to 
the detailed structure of the sound speed profile. 
Deavenport studied axial focusing using an Epstein profile, 

-2 2 

c = A sech (z/h) + B tanh(z/h) + D, (1.10) 

with solutions which are hypergeometric functions [25]. He 
found a greater degree of focusing than Williams and Horne, 
but less than from Hirsch‘s symmetric profile. 

Bucker and Morris studied sound transmission in the 
surface duct, also using an Epstein profile for the sound 
speed [26]. The results compared favorably with 
transmission characteristics derived from experiment and 
ray-tracing models. 

In 1970 Williams published an outline of a normal 
mode finite difference technique [27] which started from 
either one or both boundaries. Using a finite difference 
equation based upon the vertical Helmholtz equation, the 
method iterated across the water column to either the other 
boundary or an internal depth. An iterative approach was 
used to find the mode eigenvalues by satisfying either the 
opposite boundary condition, or continuity of acoustic 
pressure and vertical velocity at the internal depth. A 
number of other models used this technique, - notably the 
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model of Newman and Ingenito [28], and the model by Kanabis 
£29]. The Newman and Ingenito model started with a fluid 
boundary condition at the bottom and used a finite 
difference technique to compute the eigenfunction across the 
water column to the surface. The eigenvalues were searched 
by finding those values of the horizontal wave number which 
most nearly produced an eigenfunction of zero at the 
surface. The Kanabis model was similar except that it used 
a higher order finite difference equation and employed the 
adiabatic assumption for horizontal variation. 

3. Current Models 



In May 1973 the Acoustic Environmental Support 
Detatchment (AESD) of the Office of Naval Research held a 
workshop on non-ray-tracing acoustic propagation techniques 
[30]. The results from seventeen acoustic models were 
compared for calculations using a variety of sound speed 
profiles. A synopsis of the models and their results is 
given by Ref. [30]. A distinction was made between purely 
normal mode (in the sense of a complete set of orthonormal 
functions) and residue series models (where the modes are 
derived from the poles of the wave function) . At the AESD 
workshop a new model which represents a distinct departure 
from previous methods was introduced. The parabolic 
equation method, outlined by Fock for electro-magnetic 
propagation in the atmosphere [31], assumes the wave 
function to be a space dependent function, V, multiplied by 
a Hankel function. The assumed solution is applied to the 
Helmholtz equation, and terms in the second derivative of 7 
with respect to range (r) are dropped (this assumes that the 
waves are predominantly radial - that is the corresponding 
rays are at relatively shallow angles) . The following 
parabolic differential equation results: 



idV + J_d£V + JsX» (r , z) - 1 ]V = 0 , 
dr 2k Z17 2 



( 1 . 11 ) 
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where k is an average wave number and N is an index of 
refraction. Since this is a parabolic differential 
equation, with a given starting solution, and surface and 
bottom boundary conditions, the solution can be solved in 
successive range steps. The method allows full range 
dependence for the sound speed profile. 



B. OPERATIONAL MODELS 

The 0. S. Navy has a need to calculate low-frequency 
propagation loss over long ranges (in excess of 1000 km) . 
The Navy currently has variations of four standard 

operational transmission loss models and a number of 
experimental models. The operational models include the 
following: 

• FACT - The Fast Asymptotic Coherent Transmission 
Model uses modified ray theory and is primarily for 



passive sensor 


appli 


cations 


at 


low 


to intermediate 


frequencies 


(50 to 


2000 


Hz) . 


The 


FACT 


model 


computes 


propagation 


loss 


out 


to ranges 


representing 


several 


convergence 


zones 


[32]. 


FACT 


has 


been 


extended 


in an 


experimental 


version 


with 


multiple 


profiles 


and a 



bathymetry composed of linear segments along the line of 
bearing. The ray invariants of D. M. Milder [20] are 
used to connect rays between profiles. 

» NISSIM - The Navy Interim Standard Surface Ship 
Model is a ray theory model which is used primarily for 
short ranges and active sensors [33,34]. It has had some 
limited application to lower frequency passive sensors. 

■ HP-70 - This model is a multiple profile ray-tracing 
model which uses linear sound speed gradients in both the 
vertical and horizontal directions. It is used for 
passive sensors and great ranges. 

u APE - The Acoustic Parabolic Equation model, based 
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on the parabolic aquation method (Eg. (1.10)), has 
recently become an operational model. It is primarily 
used for very low frequencies and highly range dependent 
environments. 



C. SCOPE OF THIS WORK 

1 . Objecti ves 

The objective of this work is the creation of a 
computer model to estimate the transmission loss of low 
frequency CW sound over long ranges and through varying 
sound speed profiles. The model is intended to be 

relatively fast, sacrificing exactness (in the mathematical 
sense) for computational speed. The underlying theory is 
based heavily upon the WKB (Wentzel-Kramer-Brillioun) 
approximation with the eigenvalues, eigenfunctions, and 
normalization all based upon the WKB formulae. The model is 
intended for the following tasks: 

• Calculation of the detailed (fully coherent) 
transmission loss for ranges up to a few hundred nautical 
miles. In this case any sound speed variation with range 
must be slow enough to permit use of the adiabatic 
assumption. 

■ Calculation of the incoherent or averaged 
transmission loss for long ranges (thousands of nautical 
miles) . 

2. Departure fro m Previous Wo rks 

The following are a number of the major differences 
between this work and previous normal mode models: 

■ The sound speed profile is assumed to consist of 
layers (up to 50) in which the sound speed varies linearly 
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2 

in c (Eg. (1.6)). Previous models have used the analytic 
solution for such layers to form the eigenfunctions 
(modes) and find the eigenvalues (Tolstoy and May [22], 
Gordon and Petersen [30], and Bartberger and Ackler [35]). 

2 

In this model the c linear profile is primarily used for 

the ease with which the vertical wave number may be 
integrated - an operation required for many of the WKB 
calculations. The analytic solution (the Airy functions) 
is used as a mode solution only where the WKB 
approximation breaks down. 

■ The sound speed profile is divided into separate 
ducts corresponding to different families of modes. The 
calculation of eigenvalues and modes is carried out family 
by family. 

• Within each family the behavior of the eigenvalue as 
a function of mode number is approximated by a fifth-order 
polynomial. This polynomial provides the first estimate 
of the eigenvalues for each mode. For calculation of the 
averaged or incoherent transmission loss at great range, 
this estimate is sufficiently accurate for computation of 
the eigenfunctions. For detailed transmission loss at 
closer range, the eigenvalue is refined by a maximum of 
five iterations. 

■ The WKB approximation has been extended to 
incorporate the effects of a . barrier between multiple 
sound channels, and the effects of the surface and bottom 
boundaries. The eigenfunction is found by a combination 
of WKB, Airy and extrapolated functions. 

■ A range dependent environment has been included by 
the use of the adiabatic assumption. A previous 
difficulty with the adiabatic assumption for normal mode 
use has been the ambiguities for modes in different sound 
channels of multiple channel profiles. This difficulty 
has been overcome with the use of a formalism which 
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matches modes within equivalent channels and approximates 
the transition of modes between channels. 



D. REFERENCES 



1 Ewing, Maurice and J. Lamar Worzel, "Long Range Sound 

Transmission", in Pro pag ation of S ound in the Ocean , 
Geol,. Soc. Am. Memoir 27, j 5 Ocf 71348. “ 

2 Herring, C., "Transmission of Explosive Sound in the 

Sea", m Physics of Sound in the Sea, National Defense 

Research Council, Summary TecEnical Report of Division 
6, Vol. 8, 1946. Reprinted by Headquarters Naval 

Material Command, Washington, D.C., 1969. 

3 Pekeris, C. L., "Theory of Propagation of Sound in a 

Half-Space of Variable Velocity under Conditions of 
Formation of a Shadow Zone”, J. Acoust. Soc. Am., Vol. 
18, No. 2, Oct 1946, pp. 295-315. 

4 Pekeris, C. L., "Theory of Propagation of Explosive 

Sound in Shallow Water", in ProjDaaation of Sound in th e 
Ocean , Geol. Soc. Am. Memoir”27, 75“5cf 7748. - 

5 Tolstoy, Ivan, "Note on the Propagation of Normal Modes 
in Inhomogeneous Media", J. Acoust. Soc. Am., Vol. 27, 
p.274. Mar 1955. 

6 Bremmer, H. , "The W.K.B. Approximation as the First 
Term of a Geometric-Optical Series", Commun. Pure 
Appl. Math., Vol. 4, 1951, pp. 105-115. 

8 Tolstoy, I. , "Dispersion and Simple Harmonic Point 
Sources in Wave Ducts", J. Acoust. Soc. Am., Vol. 27, 
p. 897, Sep 1955. 

7 Ewing, W. Maurice, Wenceslas S. Jardetzky and Frank 
Press, Elast ic Waves in La yered Media , McGraw-Hill, 

9 Tolstoy, I. , "Shallow Water Propagation Tables, Three 
Layer Waveguides", Hudson Labs, Columbia Univ. , Tech 
Rept 75, Dec 1960, AD-256-699. 

10 Clay, C. S. , "Propagation of Band Limited Noise in a 
Layered Waveguide", J. Acoust. Soc. Am., Vol. 31, 
p.l473, Nov. 1959. 

11 Tolstoy, I.. "Guided Waves in a Fluid with a 

Continuously Variable Velocity Overlying an Elastic 
Solid: Theory and Experiment", J. Acoust. Soc. Am., 

Vol. 32, p. 81 , Jan I960; 

12 Williams, A. 0., Jr., "Some Effects of Velocity 
Structure on Low-Frequency Propagation in Shallow 
Water", J. Acoust. Soc. Am., Vol. 32, p.363. Mar 1960. 

13 Williams, A. 0., Jr., and R. K. Eby, "Acoustic 
Attenuation in a Liquid Laver over a ’Slow* 
Viscoelastic Solid", J. Acoust. Soc. Am., Vol. 34, 
p.836, June 1962. 



25 



14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 



Bucker, H. P. , "Normal Mode Sound Propagation in 
Shallow Water". J. Acoust. Soc. Am., Vol. 36, p.251, 
Feb 1964. 



Bucker, H. P., and H. E. Morris, "Normal-Mode Intensity 
Calculations for a Constant-Depth Shallow Water 
Channel", J. Acoust. Soc. Am., Vol. 38, p.1010, Dec 
1965. 



Kornhauser, E. T., and W. P. Raney, "Attenuation in 
Shallow Water Propagation Due to an Absorbing Bottom", 
J. Acoust. Soc. Am., Vol. 27, p. 689, July 1955. 



Eby, R. K., A. 0. Williams, Jr., R. P. Ryan, and Aul 
Tamarkin, "Study of Acoustic Propagation in a 
Two-Layered Model", J. Acoust. Soc. Am., Vol. 32, p. 
88, Jan 1960. 

Clay. C. S., "Effect of a Slightly Irregular Boundary 
on the Coherence of Waveguide Propagation", J. Acoust. 
Soc. Am., Vol. 36, p. 833, May 1964. 

Weston, P. E., "Guided Propagation in a Slowly Varying 
Medium", Proc. Phys. Soc. London, Vol 73, p. 365, 1958. 

Milder, D. Michael, "Ray and Wave Invariants for SOFAR 
Channel Propagation", J. Acoust. Soc. Am., Vol. 46, p. 
1259, 1969. 



Smith, P. W., "Averaged 
Range-Dependent Channels", 
55, p. 1197, June 1974. 



Sound Transmission in 
J. Acoust. Soc. Am. , Vol. 



Tolstoy, I. and J. May, "A Numerical Solution for the 
Problem of Long Range sound Propagation in Continuously 
Stratified Media, with Applications to the Deep Ocean", 
J. Acoust. Soc. Am., Vol. 32, p. 655, June 1960. 



Hirsch, Peter, "Acoustic Field of a Pulsed Source in 
the Underwater Sound Channel", J. Acoust. Soc. Am., 
Vol. 38, p. 1018, Dec 1965. 

Williams, A. 0., Jr., and William Horne, "Axial 
Focusing of Sound in the SOFAR Channel", J. Acoust. 
Soc. Am., Vol. 41, p. 189, June 1967. 



Deavenport, R. L. , "A Normal Mode Theory of an 
Underwater Acoustic Duct by Means of Green's Function", 
Radio Science, Vol. 1, p. 709, 1966. 



Bucker, H. P., and H. E. Morris, "Epstein Normal-Mode 
Model of a Surface Duct", J. Acoust. Soc. Am., Vol. 41, 
p. 1475, June 1967. 

Williams, A. 0., Jr., "Normal Mode Methods in the 
Propagation of Underwater Sound", in Underwater 
Acoustics. ed. by R. W. B. “Stephens, 
Siley-Infer science, 1970. 

Newman, A. V.. and F. Ingenito, A Nor mal Mode Computer 
Pro gram for Calculating Sound Pr op ag atio n in S hallo w 
Safer wifh an Irhtrary Velocity Prof ile7“haval Research 
Ia5 HemoranUum heporf 23BT, Jan T972. 

Kanabis, W. G., Computer Programs to Calc ulate Norma l 
Mode Propagatio n an5~A ppITc atl ons~lo the Ana lysis of 
Explosive Sounh in the BIF1 hancre,“haval Underwater 
Syst ems^enfer TechhicaI”Reporf h3T97 6 Nov 1972. 



26 



30 



Spofford, C. W. , A Synopsis of the A ESP Wo rkshop on 
Acoustic P ropa gation Modell in g Using ~Hoh- Ra y- Traci ng 
Techniques# Acoustic EnvironmentaT - Uupport Detatchment 
^SUSD) Technical Note TN-73-05, Nov 1973. 

31 Fock, V. A- , Electro ma gnetic Diffraction an 
Propagation Problem s# Pergamon Press, 1965, pp.2T3-23tf 

32 Spofford, C. i.. The FACT Mgdgl. Vol. I, AESD # Maury 
Center for Ocean Science, HTC Report 109, Nov 1974. 

33 DELETED. 



34 Watson, W. H. . and R. W. McGirr, An Active Sona r 
Perfor ma nce Prediction Model, Naval Undersea "Research 
and Ue veTopment - Center IRUC) Technical Publication 
NUC-TT-286, April 1972. 

35 Eartberger, C . and L. Ackler, Nor mal Mode Solutions and 
Co mpute r Progr ams for Un de rwa ter Sound Pro pagati on, 
Part T, RavalT" Air - Development Center Report no. 
NADC-72001-AE, 4 April 1973. 



27 



• Ip- 



II. THJSOSYi HOMOGENEOUS WAVE E QUAT ION 

la what follows it is assumed that the ocean is a 
horizontally stratified medium of constant density, with a 
perfectly flat and pressure release surface, and a flat 
bottom. The input data for most ocean acoustic models 
represent a linearly-segmented sound speed profile. The 
linear segments do not represent the detailed structure of 
the ocean sound speed profile; they are a convenience for 
interpolating between the data. In ray trace models this 
linear gradient takes on added significance because of the 
simplified ray path thus formed. Because of the analytic 
solutions available, whea using normal mode techniques it is 
convenient to approximate the sound speed structure between 

-2 

data points as linear in c . The profile will be assumed 
to consist of layers within which the sound speed varies 

- 1/2 

c(z) = c [1 - a (z-z ) ] , z < z < z , (2.1) 

i i i i+1 



where c and z are the sound speed and depth at the top of 
i i 

th 

the i layer (Fig 2) . If to is the angular frequency of the 

acoustic disturbance, then the sound speed profile can be 
represented in terms of the wave number, k. For the profile 
2 

of Eq. (2.1) k varies linearly with depth, and the gradient 
2 

of k , y , is defined by 



2 2 2 2 

k = (J /c = k - y (z-z ) , z < z < z . (2.2) 

i i i i+1 
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FIGURE 2 - Sound Speed Profile 
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The horizontal wave number k and the vertical wave number 

r 

k are related by 
z 



2 2 2 

k = k - k . (2.3) 

z r 



The vertical wave number is real for k > k , and is 

r 

imaginary for for k < k ; 

. r 



k = -i 6 , 
z 



where 



2 

6 




(2.4) 



for k < k » 
r 

The theoretical development which follows can be divided 
into five steps outlined below: 

■ The homogeneous wave equation is established, and 
from it is separated a depth-dependent equation. The 
depth-dependent equation must be solved in terms of 
eigenvalues and corresponding eigenfunctions (the normal 
modes) . 

■ Two approximate solutions to the homogeneous 
depth-dependent equation, the MKB solutions, are derived. 
One oscillatory solution corresponds to the trapped field 
in the waveguide, the other, exponentially decaying, 
represents the diffraction of sound at the edge of the 
waveguide. The WK3 solutions fail where they approach 
each other at the edge of the channel (turning points) . 
From one of the 9KB solutions a characteristic equation is 
developed, whose solutions are the modal eigenvalues. 

• A third approximate solution to the homogeneous 
equation is derived. This solution is used to connect the 
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two WKB solutions. As a consequence, for each eigenvalue 
an unnormalized eigenfunction can be determined as a 
function of depth. 

» The treatment of the homogeneous wave equation is 
completed with consideration of the effects of the surface 
and bottom boundaries. 

a The final step is the introduction of an 
inhomogeneous term (representing a point source) into the 
wave equation. The inhomogeneous equation is then solved 
by the residue series of a Green's function which is 
formed from the previously developed homogeneous 
solutions. Solution of the inhomogeneous equation allows 
determination of the mode amplitude to match a point 
source. 



A. HOMOGENEOUS WAVE EQUATION 

The homogeneous lossless wave equation for velocity 
potential V is 

V 2 ^ - 1 d£ I = 0, (2.5) 

c* Ht* 

The acoustic pressure, 
potential by 

P 

For a monofrequency acoustical disturbance with time 
dependence exp (-i ut), 

P = i up V . 

The boundary condition at the surface requires that 

4* (z=0) = 0, (2.6) 

where z is the depth below the surface. The boundary 



P, is related to the velocity 



- p d¥ 

at 
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